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Abstract 

We introduce a hybrid Cole-Hopf-Darboux transformation to relate solutions of 
nonlinear and linear second order differential equations and derive a sufficient condition 
for this correspondence. In particular we show that solutions of some nonlinear second 
order equations are related to the special functions of mathematical physics through 
this transformation. These nonlinear equations can be viewed as the "class of special 
nonlinear equations" which correspond to the linear differential equations which define 
the special functions of mathematical physics. 



1 Background 

In the context of partial differential equations the Hopf-Cole transformation [1,2] and it 
generalizations [3,4] has been used extensively to linearize some nonlinear partial differential 
equations such as the Burgers equation [5,6,7]. 

For ordinary differential equations a "similar" transformation has been used for some 
time. In fact it is well known [8] that the Ricatti equation 

ip{x)' + A(x)i){xf + B(x)tp(x) + CxO) = (1.1) 

where primes denote differentiation with respect to x, can be linearized by the transformation 

The resulting linear equation for <p(x) is 

A{x)4>{xf + (B(x)A(x) - A{x)')4>{x)' + Ci(x)A(x) 2 0(x) = (1.3) 
However the operator 

T = 4- + A(xU(x) + B(x) (1.4) 
ax 

can be used repeatedly to generate higher order differential equations that can be linearized 
by the transformation (11.21) . That is the differential equation 

T n ip(x) = Q(x) (1.5) 

where T n ip(x) = T(T"^ 1 )iIj(x) can be linearized by the transformation (jl.2p . 
In particular for n = 2 we have 

T 2 ijj = (-^- + A(x)ip(x) + B(x))(£- + A(x)ip(x) + B(x))ip(x). (1.6) 

Thus for n — 2, and v4(x) = 1, B(x) = the nonlinear differential equation (I1.5P becomes 

ijj{x)" + 3tfj(x)^(x)' + ij(x) 3 = Q(x). (1.7) 

Applying the transformation (11.21) to this equation yields 

0(x)'" = Q(x)0(x) (1.8) 
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This shows that the transformation (jl.2p relates some nonlinear equations to a linear one 
of higher order. However there is another class of transformations that are used to relate 
the solutions of two linear differential equations of the same order. These are Darboux 
transformations [9-12] (which form the basis for the well known Factorization method [13]). 
In this case the operator that relates the two equations is of the form 



Our objective in this paper is to explore the possible use of some "hybrid" form of (II. 2p 
and (II. 9j) to relate the solutions of a nonlinear equation to those of a linear one of the same 
order. In particular we are motivated by the fact that the solutions of some nonlinear 
second order equations are related to the special functions of mathematical physics through 
a transformation similar to ( 11. 21) . In a certain sense these equations define then a class of 
" special nonlinear differential equations" . Furthermore we show that Painleve II equation has 
a solution that can be expressed in terms of Airy functions for a special set of its parameters. 

The plan of the paper is as follows in Sec. 2 we present the general technique which relates 
solutions of linear and nonlinear equations. In Sec 3 we specialize to a subset of this general 
method and provide and intrinsic test for the applicability of the method. Sec 4 explores the 
relationship between some nonlinear equations and the special functions of mathematical 
physics. Sec 5 provides some additional examples of nonlinear equations whose solutions are 
related to those of a linear equation. We end up in Sec. 6 with some conclusions. 

2 A Generalized Transformation 

We shall say that the solutions of the equations 



D = C(x) + 



dx 



(1.9) 



if>(x)" = S(x) + V(x)i/i(x) + W(x)ij(x) 2 + R(x)ip(x) 3 + \i/i(x) 



(2.1) 



and 




(2.2) 



are related if we can find functions P(x) and Q(x) so that 




(2.3) 
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We observe that, in principle, the terms \ip(x), A0(x) in (12. ip . (12.21) can be absorbed by 
V(x) and U(x) respectively. Furthermore (12.11) can take the more general form 



ip{x)" = S(x) + V{x)ij){x) + VtixWx)' + W{x)i){xf + R{x)i){xf. (2.4) 

In this case we can find p(x) so that V\(x) = — 2^y. Introducing = p(x)ip(x), (12. 4p 
becomes 

C(x)» = p(x)S(x) + (V(x) + ^)C(x) + ^-m 2 + ^(x) 3 . (2.5) 

which has the same form as (12.11) . 

To classify those nonlinear equations (12. ip which can be "paired" with a linear equation of 
the form ( 12. 2p we differentiate ( 12.3P twice and in each step replace the second order derivative 
of 0(x) by U(x)4>(x) + K(x)4>(x)' + A0(x). We then use (12. ip to eliminate ip( x )"- As a result 
we find that the following equation must hold; 

a 3 (x)(f)(x)~ 3 + a 2 (x)4>(x)~ 2 + ai(x)0(x) _1 + a (x) = (2.6) 

where 

a 3 (x) = -Q(x) 2 R(x) + 2, a 2 (x) = -2Q{x)' - 3R(x)P(x)Q(x) 2 - W(x)Q(xf - 3K(x)Q(x) 

(2.7) 

ai ( x ) = Q(x)" + 2K(x)Q(x)' - (2.8) 
Q(x)[3P(x) 2 R(x) + 2U(x) + V(x) + 2P(x)W(x) - K(x)' - K(x) 2 + 3A] 

a (x) = 2(U(x) + X)Q(x)' + Q(x)U(x)' + P(x)" - W(x)P(x) 2 - V(x)P(x) - (2.9) 
R(x)P(x) 3 - XP(x) + K(x)Q(x)(U(x) + A) - S(x) 

To satisfy (12. 6p it is sufficient to let a^x) = 0, i = 0,1,2,3. We use these conditions 
to express S(x), V(x), W(x), R(x), K(x) and U(x) in terms of the parameters P(x),Q(x). 
From (12. 7p we get 



Substituting these results in (12. 8 j) we obtain an equation which we can solve for V(x) 

Q(x)Q(x)" + AP{x)Q{x)' + 6P(x) 2 - Q(x) 2 (2£/(x) - 3A) 
= + (2 ' n) 

Substituting these expressions in (12. 9p we obtain the following first order linear equation for 

U(x) 

Q(x)U(x)' + (2Q(x)' + 2P(x) + K(x)Q(x))U(x) + P{x)" + (2.12) 

2 { a - (* w + 1 gw ' + 2XP(x> - o$ (2P(x) + 3K(x)Q(x)) - 

P(x)0(x)" 

^jy\ X) + K(x)(AQ(x) - K(x)P(x)) - P(x)K(x)' - S(x) = 0. 
Q(x) 

Equation ( 12.12ft is a linear differential equation for U(x) which can be solved by standard 
methods once P(x) and Q(x) and K(x) have been specified. 



3 Solutions with Q(x) = 1 

When Q(x) = 1 we have R(x) = 2 and 

V(x) = QP(x)(P(x) + K(x))-2U(x)-3X+K(x) 2 + K(x)', W(x) = -QP(x)-3K(x). (3.1) 

Eq. ( 12.121) simplifies and we have 

U{x)' + 2P(x)(U(x) + A - P(x) 2 ) + P(x)" - S(x) + (3.2) 
(U(x) + X)K(x) - P(x)(K(x) 2 + 3K(x)P(x) + K{x)') = 0. 

Since this equation contains two unknown functions U(x) and K{x) it is natural to by split 
it into two equations 

K{x)' + K(x) 2 + (3P(x) - p^yW*) + |W = 0, (3.3) 

and 

C/(x)' + 2P(x)(f/(x) + A - P(x) 2 ) + P(x)" + K(x)f/(x) = 0. (3.4) 
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Eq. (13.31) is independent of U(x) and can be solved for K(x) once P(x) and S(x) have been 
specified. Actually (13. 3p is a Ricatti equation whose linear second order form is 

y(x)» + (3P(x) - j^)v{x)' + = (3.5) 

where K(x) = We see that by proper choice of P(x) and S(x) the solution y(x) can 

be a special functions of mathematical physics. In this case K(x) will be the logarithmic 
derivative of such functions. 

Another possible decomposition of (I3.2p is to as follows 

K(x)' + K(x) 2 + (3P(x) - -£—)K(x) = 0, (3.6) 

P[x) 

and 

U(x)' + 2P(x)(U(x) + A - P{xf) + P{xf + K{x)U(x) - S(x) = 0. (3.7) 

If we then let S(x) = — 2P(x) 3 , (13. 7p becomes linear in both U(x) and P(x). As a result one 
may choose P(x) and solve this equation for U(x) or choose U(x) and solve for P(x). 

Another option to find solutions to (I3.2p is to cancel the nonlinear terms in P(x) by 
making the ansatz 

S( x ) = -3P(x) 2 K(x) - 2P(x) 3 (3.8) 

the equation then becomes 

P{x)" + (2U(x) + 2A - K(x)' - K(x) 2 )P(x) + K(x)(U(x) + A) + U(x)' = 0. (3.9) 

which is a linear equation for P(x) once U(x) and K(x) were chosen. 

In the following we use all these strategies to to find nonlinear differential equations 
whose solutions are related to those of a linear differential equation by the transformation 

We now consider the following practical question: Suppose one is considering a nonlinear 
differential equation of the form given by (12. ip . Under what conditions one can find a linear 
differential equation (12. 2p whose solutions are related to it by (12. 3p with Q(x) = 1. 

Theorem: A sufficient condition for the solution of (12. ip to be related to an equation of 
the form (E2D by (E3D with Q(x) = 1 is that R(x) = 2 and 

S(x) = -\ -V(x)' + - [-W(x)" + (V(x) + W(x)' + \)W(x)\ - \w{xf \ (3.10) 
2 1 3 27 j 



When this condition is satisfied one can choose K(x) = and 

U{x) = - l -V{x) + ^W{xf-\\ (3.11) 

We observe that the condition ( I3.10P is an intrinsic condition on the coefficients of ( 12. ip . 
Proof: From (13. ip we have 

P(x) = -\w(x) - \k(x) (3.12) 
o 2 

Substituting this expression in the formula for V(x) it follows that 

V(x) - \w(xf + \k(x) 2 + 2U(x) + 3A - K(x)' = (3.13) 
6 2 

Solving ( I3.13P for U(x) we find that 

U(x) = -~V(x) + ±W(x) 2 - l -K{xf - + l -K{x)' (3.14) 

Finally substituting these expressions for P(x) and U(x) in (13. 2p the condition ( I3.10P follows. 
If we let K{x) = in ( I3.12p . (l3.13p and ( I3.14p the expression for U(x) reduces to the one 
given by ( 13. lip . ■. 



4 Relationships to the Special Functions 

In this section we explore the relationship between some nonlinear equations and the function 
of mathematical physics. The purpose of our treatment is to highlight this new relationship 
and is not comprehensive. In all cases we let Q(x) = 1. 

Case 1: Exponential and Trigonometric Functions: 

We start by considering the harmonic oscillator equation, 

(f)(x)" = -cu 2 (p(x), co^O. (4.1) 
For this equation U(x) = —u> 2 , K(x) = 0, and A = 0. Hence from ( 13 .9p we have 

P{x)" - 2oo 2 P(x) = 0. 

Therefore 

P( x ) = C x e^ x + C 2 e-^ x 
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Computing V(x) W(x) and S(x) using (13. ip and (13.81) we finally obtain the following differ- 
ential equation for ip{x) 



ij)(xf - (6P(x) 2 + 2u 2 )ip + 6P(x)^(x) 2 - 2ip{x) 3 = 2P(xf (4.2) 

whose solutions of are related to those of (14. 2p by (12. 3p . 
If we consider 

(j){xf = u 2 <p(x) = 0, co ^ 0, (4.3) 

then 

P(x) = C\ sin(\f2ux) + C 2 cos(\/2a;x) 
and the differential equation for ip(x) becomes 

ij){x)" - (6P(x) 2 - 2w 2 )^(x) + QP(x)^(x) 2 - 2ifj(x) 3 = -2P{xf (4.4) 

For the case where u = i.e. 

y(x)" = (4.5) 

then 

P(x) = Ci + C 2 x. 

and the equation for ip(x) is 

ip(x)" - 6P{x) 2 ip{x) + 6P(x)ij(x) 2 - 2ip(x) 3 = -2P(xf (4.6) 

Similar treatment can be made for the general second order equation with constant 
coefficients. We omit the details. 

Case 2: Legendre polynomials 

The differential equation for the Legendre polynomials is 



0(*y = -^<K*)' - ^±i)#s). (4.7) 

1 — x z 1 — x z 



That is U(x) = — vjzd^l ; K{x) = and A = 0. Substituting this in (I3.9P we find that 
particular solution for P(x) is 

/ x 2n(n + l)x 

P(x) = -— V ./ -, n ^ 1 4.8 

v ; (n 2 + n - 2)(1 - x 2 ) r V ; 
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a 



(The solution for n = 1 is available but we shall not elaborate on it further). Using (13. ip we 
have 

A8n(n + l)x 2 2n(n + 1) 2(2x 2 + 1) 

K ) (n + 2) 2 (n- 1) 2 (1 -x 2 ) 2 1-x 2 (1 - x 2 ^ 2 ' 1 
From (13. 8 p and (13. ip we obtain 



8n> + l) 2 (n + 3)(n-2)x 3 
W (n + 2) 3 (72-l) 3 (a; 2 -l) 3 ' 1 J 

= -QP(x) - 3K(x). (4.11) 

For these functions S(x), V(x), W(x) (and R(x) = 2) the solutions of (12. ip are related to 
the solutions of Legendre equation by the transformation (12. 3p . 

Case 3: Bessel Functions 

The differential equation for Bessel functions is 

m* = - (4.i2) 

Jb Jb 



That is U(x) = — X 2 P , K(x) = — - and A = 0. Substituting these in (13. 9p we find that the 
general solutions for P(x) with p — 0,1 respectively are 

z=v - 2Xt (413) 

x 2x 

w=w+ ^.j«ai)^M + i (4 , 4) 

2x 2x 



oo 



where 

r(0,z) = / r^-'tft 
If we choose for n = the special solution PoO^) — ^> we find that 

The differential equation for ip(x) is 

v^r = ^3 + ( 2 + 2^ + 2 ^ x ) 3 - ( 4 - 15 ) 

The solutions of this nonlinear equation are related to the Bessel functions of order zero by 
the transformation (12. 3p . 
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Case 4: Hermite polynomials 

The differential equation for Hermite polynomials is 

(f)(x)" = 2x(p(x)' - 2n<f>(x). (4.16) 

Hence U(x) = —2n, K(x) = 2x. Substituting these in (13.91) we find that an explicit particular 
(and general) solution for P(x) in terms of elementary functions can be obtained for even n. 
For n = 2, 4 these particular solutions are 

8 [tt(1/4 + x 2 ) exp(x 2 )er/(x) + ^x] 

^2{x) -j= , l 4 - 1 '; 

'7T 



64 [tt(3/16 + 3/2x 2 + x 4 ) exp(x 2 )erf(x) + y¥x(l + x 2 )] 
r^4{X) — , (4.1oJ 

3V 71 " 

where 

/ exp(-t 2 )rft 
er/(x) = 2^- 



'7T 

The computation of the functions /S'(x), V(a;), W(a;) using ( 13. ip and ( 13. 8p is straightforward. 
Case 5: Laguerre Polynomials: 

The differential equation for the Laguerre polynomials is 

1 — T 11 

Mx)" = -^-t^ x y - r l^ x ). (4.19) 

x x 

That is U(x) = —-, K{x) = — and A = 0. Substituting these in ( 13. 9 p we find that a 
particular solution for P(x) with n = 2 is 

PJ X ) = -r(0, -x)e~ x x + e s r(0, + 2x - 2) - 2 + - (4.20) 

X X 

We omit the computation of V(x), W(x) and S(x) using ( 13. ip and ( 13.81) which is straight- 
forward. 

Case 6: Painleve II Equation 

Painleve II equation is 

ip(x)" = 2^(x) 3 + xtp(x) + a (4.21) 

where a is a constant. In the notation of the previous sections V(x) = x, W(x) = 0, 
R(x) = 2, S(x) = a and A = 0. With Q(x) = 1 (14T20|) satisfies the constraint (13401) when 
a=-|. 
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With K(x) = 0, (13. lip and (13. 12)) yield under the present settings P(x) = and 

U(x) = -\v{x) = -\ (4.22) 
i.e the differential equation for 4>(x) is 

<j>{x)" = ~<j>{x) (4.23) 

whose general solution is 

4>{z) = C X M{~) + C 2 Bi{-^) (4.24) 

where Ai(x), Bi(x) are the Airy wave functions. This solution is related to the solution of 
( 14 .21)) by the transformation ( 12 .3)) (with P(x) = 0). We observe that other solutions of this 
equation are possible if we let K(x) ^ 



5 Some Special Cases 

In this section we present some explicit solutions to the equations which pair the solution of 
a nonlinear equation with a linear equation using the algorithm which was presented in the 
previous sections. In all cases we let Q(x) = 1. 

Example 1: In this example we let K(x) = and S(x) = and use (1 3 . 2 p . Choosing 
P(x) = \ we obtain for n = 1 

U{x) = C l x- 2b + h ^-^-\ (5.1) 

x z 

and the corresponding expression for V(x) is 



v (l) = 2ClI - + ?(?M-, 

X z 



For n / 1 we have 

^) = ^exp^-^j + - TT , + --A, (5.2) 

and 

fb\ 2n ( 2bx^ \ 2bn 26 2 
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For n = 1 the equation for (f>{x) with C% = is 



_ ^±I1 0(X ) = o, 



(5.3) 



whose general solution is 

<f>(x) = D x x b+1 + D 2 x- b . 

The corresponding nonlinear differential equation for ip(x) is 

. , ,„ 26(26 + 1) , . . 66^(x) 2 n , . , 3 n 
^(a;)" ^5 -ip{x) H Z±J- _ 2^(x) = 0. 

In this case the explicit relationship between <p(x) and i/j(x) (as postulated in ( 12.31) ) is 

x 0(x) 

It is straightforward to verify that this is actually a solution of (15. 51) . 



(5.4) 



(5.5) 



(5.6) 



Example 2: Let P(x) 



-2a; 



S'(x) 



-4na; 



and A = 0. Eq. ( 13. 5ft becomes 



?/(x)" — 2xy(x)' + 2ny(x) = 



(5.7) 



which is the differential equation for the Hermite polynomials. Hence 

HJx)' 



K(x) 



HJx) 



(5.* 



where HJx) is Hermite polynomial of order n. Substituting these results in fl3.4[) leads to 
the following general solutions for U(x) 



U n {x) 



d-J §x 3 Hjx) exp(-^)dx 



exp(^) 



HJx) 



(5.9) 



The integral in (15. 9p can be computed explicitly for different values of n. For H (x) = 1 we 
have 



2r 2 3 
U (x) = C 1 eM—)--(3 + 2x 2 ) 



For H 2 (x) = Ax 2 - 2 



U 2 {x) 



C x exp(^) - \(Ax A + 10x 2 + 15) 
4x 2 -2 
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etc. 

From these expressions for K(x) and U(x) we find that 

i v 8x 2 , * HJx)" - AxHJxY , > 3# n (x)' 

We conclude that the solutions of the nonlinear equation ( 12. ip with R(x) = 2 and V^x), W n (x) 
in (15.101) are related to those of the linear differential equation (12. 2p with U(x) and K(x) 
given by (15. 9p and (15. 8p by the transformation (12.31) . 

Example 3: In this example we let K(x) = 0, P(x) = a + bx, S(x) = — 2P(x) 3 and use 
(13. 7p to compute U(x). This yields 

U(x) = Ci exp[-(2a + bx)x] - A, V(x) = 6(a + bx) 2 - A - 2Ci exp[-(2a + fcx)x]) (5.11) 

and W(a;) = -6P(x) (from (12~T0D ). 

Example 4: Here as in the previous example we let K(x) = 0, S(x) = — 2P(x) 3 and 
A = 0. However we now let U(x) = — -% and use (13. 7p to compute P(x). For d^lwe obtain 



For a 



P( x ) = - " + + C^ 1 / 2 "^ /3 = v/TTS^. (5.12) 

a — 1 x 



. , 6 In x + 2 „ 9 Co / % 

P(x) = — + Cix 2 + — . 5.13 

X X 



The corresponding V(x) in this case (viz. a = 1) is 

V(x) = -^ + P(x) 2 . (5.14) 

x z 

If we reverse the roles of P(x) and U(x) i.e fix P(x) = — -% and use (13.71) to compute 
[/(x) we find 



[%) = CVexp(— -)-3(^- — + 2^)- (5.15) 



x \ar ax 

Example 5: Under the same settings of example 9 we specify U(x) = 2exp(2ax) and 
use (13. 7p to compute P(x). We find 

P(x) = C 1 J (z)+C 2 Y (z), z= 2eXp(ax) (5.16) 

a 

where Jq and Y are Bessel functions of order zero of the first and second kind. The corre- 
sponding values of V(x) and W(x) are 

V(x) = QP(x) 2 - Ae 2ax , W(x) = -6P(x). 
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6 Conclusions 



We demonstrated in this paper that the hybrid " Cole-Hopf-Darboux" operator ( 12. 3p can be 
used to relate the solutions of some nonlinear and linear second order differential equations. 
The algorithm is straightforward to apply and we presented several examples which demon- 
strated the scope of the method and its potential interest for problems in mathematical 
physics. From a perspective those nonlinear differential equation whose solutions can be ex- 
pressed in terms of the special functions of mathematical physics (through the transformation 
(12. 3 p should obviously be viewed as a new " special class of nonlinear equations" . 

Research into the possible application of the transformation ( 12.31) to partial differential 
equations is on going. 
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